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A thermodynamic theory of ferroelectric thin film multilayers is developed. The free energy func- 
tional is written down using a multilayer model in which c-domain layers of one ferroelectric material 
alternate with a-domain layers of a second ferroelectric material. We assume that the interfaces are 
perfectly sharp and that the polarization at these boundaries is zero. The renormalization of the 
free energy coefficients due to the stresses in the films and to the depolarizing field was taken into 
account, as well as the renormalization of the coefficients of the polarization gradients. The equilib- 
rium inhomogeneous polarization temperature and its thickness dependence were determined from 
the solutions of the Euler-Lagrange equations resulting from the minimization of the free energy 
functional. A thickness induced ferroelectric phase transition is shown to exist and its transition 
temperature and critical layer thickness depend on the domain orientation. The criteria for "a/c", 
"c/c" and "a/a" domain multilayering are calculated and expressed via coefficients of the free energy 
density and the layer thickness. 

The differential equation to calculate the inhomogeneous dielectric susceptibility was found to be a 
Lame equation. The calculated multilayer susceptibility diverges at the transition temperature of 
the thickness induced ferroelectric phase. This divergence is shown to be the origin of the giant di- 
electric response observed in some multilayers. The theory gives an excellent fit to the temperature 
dependence of the giant susceptibility observed recently in multilayers of PbTiOs-Pbi-jjLajjTiOs 
{x = 0.28). The equation determining the time and space dependence of the polarization is pro- 
posed and solved. The calculated dispersion law for the nonlinear polarization waves in multilayer 
structures reveals a critical wave vector for a given layer thickness, or a critical layer thickness for 
a given wave vector, for which the frequency u = 0. The frequency of the polarization waves is also 
shown to increase with layer thickness. 



PACS numbers; 6865.+g, 7755.-ff, 7784.-S 
I. INTRODUCTION 

Artificial thin film multilayers composed of alternat- 
ing layers of different materials have been the subject 
of study for many years. Much attention has been de- 
voted to semiconducting metallic magnetic ||] 
and superconducting 0| multilayer structures. Increas- 
ingly frequent investigations of ferroelectric oxide mul- 
tilayers are now taking place due to their technologi- 
cal promise: the artificial modification of structural and 
physical properties for use in dielectric capacitors, mem- 
ory systems, pyroelectric detection and other types of 
devices makes them technologically attractive (see e.g. 
1^,^ and references therein), and there is also a ba- 
sic interest in creating model structures for which to 
study fundamental questions related to ferroelectricity 
on increasingly smaller length scales. Recent experi- 
ments have revealed many unusual properties, such as gi- 
ant dielectric response (tI , anomalies in the ferroelectric- 



paraelectric phase transition temperature |^ , in superlat- 
tice growth and in superlattice structural anomalies 
10 . In the latter work, investigation of PbTiOa /BaTiOs 
(PT/BT) multilayers with bilayer PT/BT thickness be- 
tween 50 A° < A < 360^1° 10 shows that "c" -domain BT 
layers alternate with "a" -domain PT layers. This is in 
contrast both to the "c" -domain structure exhibited by 
thin single PT films of the same thickness grown on the 
same substrate. Raman spectra measurements reinforce 
this x-ray determined orientational anomaly in addition 
to revealing the existence of a superlattice wavelength 
dependent mode, whose frequency increases with increas- 
ing A. Up to now there have been few theoretical studies 
of ferroelectric multilayer structures. Calculation of the 
dielectric response has shown that the spatial distri- 
bution of the layer thickness can lead to enhanced dielec- 
tric properties over a broad temperature range. A phe- 
nomenological theory has recently been developed for 
thin film multilayers. The spontaneous polarization and 
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the dielectric susceptibility were calculated numerically 
with parameters appropriate to PT/BT. They showed 
that a thickness induced ferroelectric phase transition oc- 
curs only when the strength of the interfacial coupling is 
weak. 

In this work we will consider a multilayer built up from 
alternating (100) and (GDI) ferroelectric thin films (such 
as PbTiOs and BaTiOa) epitaxially grown in the cubic 
paraelectric phase onto a cubic (001) substrate. We use a 
thermodynamic approach to describe the static and dy- 
namic multilayer properties. The static properties are 
considered in section II in which section IIA is devoted 
to the description of the free energy of a multilayer with 
contributions coming from the internal stresses and the 
depolarizing field. Calculations of the inhomogeneous po- 
larization and the thickness induced ferroelectric phase 
transition are performed in IIB. The criterion for the 
presence of "a/c", "a/a" or "c/c" layer structures in the 
superlattice is determined in IIC. The differential equa- 
tion giving the static dielectric susceptibility of multi- 
layers is considered in IID. The dynamic properties are 
treated in Section III. In IIIA the differential equation 
governing the time and space dependence of the polar- 
ization is written down and solved, and in IIIB the dis- 
persion law for nonlinear waves formed by the polariza- 
tion in the multilayer is reported. In the Discussion we 
show that the theory is in good agreement with the ob- 
served giant dielectric response in PLT/PT ferroelectric 
thin film multilayers and in the Conclusion we discuss 
the applicability of our results and their future develop- 
ment. The Appendices contain the detailed calculations 
of (1) the criteria for "a/c", "a/a" and "c/c" layer struc- 
ture of the superlattice, (2) the solution of the differential 
equation for the inhomogeous dielectric susceptibility, (3) 
the multilayer susceptibility, (4) the value of the suscep- 
tibility for the thick film ferroelectric phase and for the 
thin film paraelectric phase, and (5) the Curie- Weiss law 
in the vicinity of the thickness induced ferroelectric phase 
transition. 



II. STATIC PROPERTIES 

A. Free energy 

We consider a multilayer (see Figurel) built up from 
"A" and "B" ferroelectric layers extending from z = 
—L/2 to L/2. Each layer has thickness U (i = A,B) 
so that the total multilayer thickness is L = N{Ia + Ib) 
where A = {I a + Ib) is the multilayer wavelength and 
N is the number of wavelengths in the superlattice. The 
free energy of the multilayer system can be written as 
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{^A{z) + ^B{z) + ^AAz))dz 



(1) 



where $yiand $b are the free energy densities of the A 
layers and the B layers respectively and $a,_b is the free 
energy density resulting from the interaction between the 
layers. 

We will start with the following forms for $^ and $ b , 
which can be obtained from those free energies having 
cubic symmetry, and which allows for the symmetry low- 
ering that is related to the nonequivalence of the z and 
the x, y polarization components in the film due to the 
contribution of the mechanical strains and to the depolar- 
izing field. These effects are incorporated into the renor- 
malized coefficients. 
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Since we wish to take into account the change of Px^y 
along the z direction, we have incorporated the additional 
gradient terms dPx.y/dz into Eqs.(2). 

Equations (2a) and (2b) correspond to the general case 
in which the x,y, and z polarization components exist 
in both the A and in the B layers. Experimental data 
for multilayers usually show that c-domain layers alter- 
nate with c-domain layers |^], but recent experimental 
results on PbTi03/BaTi03 multilayers have shown |ic| ] 
that c-domain BT layers alternate with a-domain PT 
layers in the multilayer structure. In order to be able 
to treat this interesting latter case we will assume that 
the A layers consist only of c-domains (i.e. Pa — Paz, 
Pax = PAy = 0) while the B layers consist only of a- 
domains (i.e. Pbx = PBy 7^ 0, Pbz — 0). The variation 
of the polarization in the A and B layers in the vicinity of 
the 2 = boundary is schematically depicted in Figure 
2 in which the natural boundary condition for this con- 
figuration is Paz{z = 0) = Pbx{z = 0) = . Since Fig. 
2 represents one modulation period of the superlattice, 
this same boundary condition must be valid for all the 
interfaces that make up the multilayer structure, namely 
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PAz{zj) = Pbx{Zj) = 0, 



(3) 



j = 0, l,2...iV- 1. 
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The superlattice periodicity also implies the following pe- 
riodic polarization condition, 



PAziz+jilA + ls)) ^Paz{z) 

Pbx{z + j{lA + Ib)) = Pbx{z) 



(4) 



Note that even in the case of a "c-domain/c-domain" or 
a "a-domain/a-domain" multilayer structure, zero polar- 
ization at the boundaries would be the most probable 
since the situation at the interfaces is physically similar 
to that of domain boundaries in bulk ferroelectrics. 

We will, in this work, neglect the interaction energy 
term ^ab due to the polarization interaction between 
the layers since this interaction decays rapidly, and be- 
cause of the boundary conditions, we will also neglect 
any surface energy contributions. Hence we can rewrite 
the free energy density (2) equations as 
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with the corresponding boundary conditions for P, a, P' , 
a' . Note that Eqs. (7a) and (7b) can be written for all 
components of Pi and aij respectively. 

In what follows we will calculate P{z) on the basis of 
Eq.(7a) and a from first term in Eq.(7b) because we will 
here make the simplifying assumption that the interfa- 
cial stresses are homogeneous and so we will neglect the 
second term in Eq. (7b). For films, the coefhcients of Eq. 
(2) are renormalized coefficients. The stresses as well as 
the depolarizing field act on the coefficients of the free en- 
ergy of cubic symmetry. In general the thickness of the 
substrate is much larger than that of the multilayer struc- 
ture. Taking this into consideration for the equilibrium 
condition of the mechanical forces, it can be shown 
that the stresses induced by upper layers in the underly- 
ing films are negligibly small. Thus the nonzero stresses 
are axx = <^yy and a^y = while azz — <^xz = '^yz ~ be- 
cause of the existence of the unstressed free surface on top 
of the multilayer. The renormalization of the free energy 
coefficients a^, a.y, (z, j, k - Voigt notation) of bulk cu- 
bic symmetry by these stresses was previously performed 
p^ , p^ for a single thin film on a cubic substrate. Us- 
ing the results of jl6| along with the depolarizing field 
contribution, the coefhcients as and a\ in Eq.(2a) can be 
written as : 
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In Eqs. (5) we have conserved what is most essential for 
the analysis of thin film multilayers: the gradient terms 
which take into account the polarization change along 
the z direction (see Fig. 2). Thus, in our model, the free 
energy density f{z) = fA{z) + fsiz) will replace ^a{z) + 
^b{z) + ^ab{z)- The most important feature of the free 
energy density of thin films is its dependence on z. In 
the most general case both the polarization P and the 
mechanical stresses a can be inhomogeneous so that their 
gradients P' and cr'should also be taken into account. 
Because of this, the free energy $ is a functional of P, 
P', a, a' , which we will write as: 
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where, for the sake of clarity, we have omitted the vector 
and tensor component notation. 

The equilibrium values of polarization and stresses or 
strains must then satisfy the Euler-Lagrange equations 



a/ 

dP 



dz dP' 



(7a) 



ai = a — X 



Qn + Qi2 

^11 + ^12 



(8b) 



Here a = aQ{T ~ T^), T^, Qfj and S(- are respec- 
tively the coefficient of the free energy, the ferroelectric 
phase transition temperature, the electrostriction con- 
stants and the elastic modulus of the A material. The 
strain X — Xi — X2 can be represented as 
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where Xmf = (6 - a)/a, Xth = (as ~ ola){T ~ Tg) and 
Xdis are respectively the misfit, thermal and disorder in- 
duced strains. In these expressions, 6 and a are the cubic 
lattice constants, ua and as are the thermal expansion 
coefficients of the two materials in the cubic bulk phase 
and Tg is the growth temperature. We point out that 
the relaxation processes related to misfit dislocations, do- 
main structure appearance and impurity diffusion pro- 
cesses can decrease X, but full relaxation can only be 
achieved in bulk materials |jT^ . The last term in Eq. (8a) 
originates from the depolarizing field Ed whose contribu- 
tion equals to —l/2EdPAz with Ed — —^ttPaz/sa for a 
free standing film without electrodes {ea is the dielectric 
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constant of the A material jT^]). The renormahzed coef- 
ficients in Eq.(2b) can be simply obtained from Eqs.(8) 
by the substitution as — s- 63, ai — > bi and the index A — s- 
B. 

The renormalization of the transition temperature in 
the layer follows from Eqs.(8) and its analog for the B 
layer. For the A and B layers the transition temperatures 
can be written as: 
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where Tcz and Tex are respectively the transition tem- 
peratures for the appearance of ^ (c-domain struc- 
ture) and Px ^ Q (a-domain structure). It is seen that 
the depolarizing field contribution decreases the transi- 
tion temperature of the c-domain layer whereas the influ- 
ence of strain on this temperature depends on the signs 
of the electrostriction constants and whether the strains 
are tensile {X > 0) or compressive {X < 0) . 

We have omitted terms to the sixth power and above in 
the polarization, which is the correct procedure only for 
second order ferroelectric phase transitions. In general, 
bulk ferroelectric materials such as BaTiOa and PbTiOa 
undergo first order phase transitions, but calculations on 
BaTiOa and PbTiOa thin films have shown ^ that the 
renormalized coefficients of and have positive val- 
ues and thus the ferroelectric thin film transition is sec- 
ond order rather than first order. For this reason we 
consider the free energy forms (2) and (5) to be those 
appropriate for multilayer films. 

To the best of our knowledge there has been no con- 
sideration of how the stresses will renormalize the coef- 
ficients of the gradient terms in the free energy of cubic 
symmetry. We have performed calculations similar to 
those in for which we have added the additional 

term AF to the cubic symmetry free energy contained in 
Ref.l6: 
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where cti = a^x — o'yy 7^ is the nonzero homogeneous 
stress in the film. Minimization of the free energy with 
respect to ai gives the renormalized coefficients: 
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Obviously, the same relations hold for the free energy 
and the (3 coefficients of the B material. The parame- 
ters 7, Sij, 6ijk are material parameters. The coefficient 
a44 is proportional to X, i.e. 0:44 7^ in films having 
nonzero homogeneous strain X (8c). The coefficients 6ijk 
are components of a sixth rank tensor which is the lowest 
rank (and so having the largest components) which will 
relate the squared gradient terms and the components of 
the stress tensor. 



B. Polarization 

The Euler-Lagrange equation (7a) with f ^ Ja + /b 
(see Eqs.(5)) for P = Paz or Pbx and P' = {dPAz)/{dz) 
or [dPsx) / {dz) makes it possible to find the equilibrium 
values of the polarization on the basis of the equations: 
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The above equations have to be solved subject to the 
periodicity conditions (3) and (4). We shall demonstrate 
the solution for Paz only, since Pbx can be obtained from 
Paz by substituting the corresponding coefficients. To 
integrate Eq.(12a) we let {dPAz)/{dz) = g{PAz)- This 
gives 



d'^PAz , dg{PAz) 

Substitution of Eq.(13) into (12a) leads to 

MPaz) 
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which gives after integration 

o-sPaz + o-ssPaz = assg^iPAz) + C33 



(13) 



(14) 



(15) 



To obtain the constant C33 we introduce the maximum 
polarization in the layer PAzm which satisfies the condi- 



tion 



dz 



|p^^=p^^„= 0, and we find 

C33 = «3-Pl2m + O'SsPazt 



(16) 



Substitution of Eq.(16) into Eq.(15) leads to 
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We now introduce the following parametrization 

Paz{z) = PAzm sine Aiz) (18) 



The relation between PAzm and Iaz can be found from 
the periodicity condition of the elliptic sine function (see 
MM), namely 
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where P^^q = — ^3/(2033) is the homogeneous polariza- 
tion in a thick film - when the derivative in Eq.(12a) 
can be neglected. Note that this is not the polarization 
in the bulk material because the parameters 03 and 033 
have been renormalized by the stresses in the layers (see 
Eqs.(8)). Because the homogeneous polarization corre- 
sponds to the mean field approximation, PazO can be 
considered as the film polarization calculated in this ap- 
proximation. Substituting Eqs. (18) and (19) into (17) 
gives: 
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where we have introduced the dimensionless variables 
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It follows from the theory of elliptic functions 19-2l|] 
that for the relation (21) the function Oa{z) has the form 
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Z3 + L,/2^ 2jlAz 
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Hence from Eq.(18) we obtain the solution in terms of 
the elliptic sine function sn: 



p . ^ p Z3 + LJ2 - 2jlAz , 

PAz{Z3) = PAzmSn , KAz 



(23) 



It can be seen that Paz{z) satisfies Eq.(12a) and the pe- 
riodicity conditions of Eqs. (4). Paz{z) is depicted graph- 
ically in Fig. 3. 
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where 
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is the complete elliptic integral of the first kind. 
Substitution of (24) into (23) gives 
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We plot the dependence Eq. (25) in Fig. 4. Since we can 
write the polarization ratio as PAzmf-^AzO — "^^Az/i^ + 
^Az) (^^6 Eq.(19)), Fig. 4 makes it possible to obtain the 
dependence of this ratio on the film thickness. We can see 
that for thick enough films (for the dimensionless length 
Iaz « 10 — 12) Pazui ~ PazO: and so this polarazation 
ratio (see the additional scale in Fig. 4) can be used as a 
criterion for distinguishing between thick and thin films. 
Moreover, Fig. 4 shows the existence of a critical layer 
thickness Iaz = such that the spontaneous polariza- 
tion in the layer will exist only for Iaz ^ t^- Thus our 
calculations yield a thickness induced ferroelectric phase 
transition which was previously discussed for multilayer 
films and for single ferroelectric thin films |p2| , p3| . 
The thickness induced phase transition temperature 
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where T^, the renormalized layer temperature, is given 

^/T^ as a func- 



by Eq(9a). The reduced temperature T^i/ j-co 
tion of reduced thickness Ia/Iao is plotted in Fig. 5. The 
characteristic thickness Ia = Iao at which t4 = is: 
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The Fig. 5 plot is in agreement with available experimen- 
tal data on ferroelectric KNb03 /KTa03 superlattices 
for intermediate multilayer wavelengths. The range of 
the existence of the thickness induced phase transition is 
given by 



Iao <Ia< 



(26c) 



with Ia 
temperature 



00 at 



T^, the renormalized transition 



Below a certain thickness, Ia < Iao there 
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is no thickness induced phase transition because be- 
comes negative. The thickness dependence of (see 
Eq.(26a)) means that the thinner the film, the weaker its 
ferroelectricity. 

The distribution of the polarization parallel to the mul- 
tilayer growth axis, Pbx{z) ,can be obtained from (23) by 
simply replacing the coefficients (see (12a), (12b)): 



as ^ 6i, 2033-^611, 033 ^ /333 
The solution for Pbx{z) thus has the form 
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Fig. 3 of course also represents Pbx{z)- 

The temperature of the thickness induced ferroelectric 
phase transition can similarly be obtained from (26), (27) 
with the substitution T/ 



rjiB rjiB 

cl ~ CX 
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Since the parameters in Eqs.(26) and (30) are different, 
the critical characteristics of the thickness induced phase 
transition should also be different in the A and the B 
layers. This may open up the prospect of engineering 
new multilayer materials constructed with several ferro- 
electric thin films (including superstructures consisting 
of several thin films in its unit cell) with a broad distri- 
bution of the transition temperature. This will result in 
a distribution of the material properties and also in any 
anomalous behavior, which may then be exploitable for 
device applications. 



situation is energetically favorable compared to a- only 
or c- only layers throughout the multilayer structure (a/a 
and c/c layering respectively). The requirement for the 
a/c multilayer to occur is that the free energy should be 
less than that for c/c or a/a multilayers. In supposing 
that the A layers have "c" -domain structure, the crite- 
rion for the existence of an a/c multilayer can be written 
as 



FBiPBx) < Fb{Pbz) 



(31) 



where Fb can be obtained by integration of the free en- 
ergy density (5) over dz and summing over the layers. 
When the inequality (31) is not satisfied we will have 
the "c/c" domain criterion. Eq.(31) will correspond to 
"a/a" domain criterion if we suppose that the A layers 
are a-domain. To calculate Fb{Pbx) and Fb{Pbz) we 
have to substitute Eq.(23) into Eq.(5a) (with A ^ B) 
and Eq.(28a) into Eq.(5b) and then perform the inte- 
gration over dz. This integration and summation (see 
Appendix 1 for details) yield the following criterion for 
"a/c" domain structure: 
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is the complete elliptic integral of the second kind. The 
function tf{k) is calculated numerically and the results 
are shown in Figure 6. One can see that ^p{ki) slowly in- 
creases for 0.3 < ki < 0.8 and for 1 > fci > 0.8 it changes 
much faster. In the scale of dimensionless thickness Iz the 
region of slow increase corresponds to 5.5 > lz,x > 3.5. In 
the thick film limit {kBi = ^) f — 3/8 and the criterion 
(32) transforms into 



C. Criterion for "a/c", "c/c" and "a/a" domain 
structures 

Up until now we have considered a multilayer built up 
of layers in which the polarization alternates between be- 
ing in the plane of the film (a-domain layers) and perpen- 
dicular to the plane of the film ( c-domain layers). This 
choice was motivated by the experimental results of Ref 
We will now look at the conditions for which this 
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This can be rewritten as 
{Tfz-T? 



< 



'J33 



bii 



(34) 



(35) 



where the transition temperatures in the thick films ^ 
are given by Eqs.(8). It is seen that for > T/^ the 
preference is for a-domain orientation in the B layers (and 
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hence for a "c/a" domain multilayer). In comparing the 
formulas for T^z and Tex (Eqs. (9a) and (9b)), we see 
that for tensile strain {X > 0) and for Q12 < and 
Qii + Q12 > 0, which is the case for many ferroelectrics 
with the perovskite structure, Tex will be higher than Tcz- 
However for compressive strain {X < 0) it is possible for 
Tex to be smaller or larger than Tcz depending on the 
depolarizing field. T^^ > T^^ can result from the depo- 
larization field's negative contribution to Tcz- However, 
when the depolarization field can be neglected (e.g. a 
film with electrodes having high conductivity p8[|) com- 
pressive strains will lead to Tex < Tcz- 

We point out that the criterion given by Eq.(35) cor- 
responds to considering the film in the mean field ap- 
proximation with homogeneous polarization P — Pq (see 
Eq.(19)). It is the function ip which takes into account 
the contribution of inhomogeneous polarization related 
to the gradients in the free energy. Thus for the general 
case, one has to use criterion Eq.(32) for "a/c" domains 
in the multilayer structure, keeping in mind that the op- 
posite condition (sign ">" substituted for sign "<" in 
Eq.(32)) implies the existence of a "c/c" -domain mul- 
tilayer structure. The criterion for "a/a" domain coin- 
cides with Eq.(32) for the case when the A layers are 
"a" -domain. The criterion depends on the free energy 
parameters, the transition temperatures, and the ratio 
of the maximum polarization in a layer P„i to the thick 
film polarization Pq. The latter in turn depends on the 
film thickness, the coefficients of the gradient terms and 
the free energy constants 63, 61. With the help of Fig. 4 or 
using the analytical formulas Eqs. (25) and (28b), Eq.(33) 
can be rewritten in terms of a dimensionless layer thick- 
ness Ibx.z which is presented as the second scale in Fig. 6. 
By keeping in mind the analytical form of Ibx.z (Eq.29) 
one can see that the criterion (Eq.32) depends on the 
film thickness and the temperature dependent free en- 
ergy parameters. Knowledge of these parameters will 
make it possible to calculate completely the criterion for 
a multilayer domain structure for a given layer thickness. 

D. Static dielectric susceptibility 

We will now consider a multilayer system in an exter- 
nal electric field E — Ez- In so doing we must add the 
term —EzPaz to the free energy density of the ferroelec- 
tric phase. From the general form of the Euler-Lagrange 
equation (7a), it follows that we have to add —Ez to 
Eq.(12a). Therefore this equation now has the form: 

2a^PAz + 4a33PL - Ez - 2«33^^ = (36) 

The applied electric field induces an additional homoge- 
neous polarization AP^^ = XzzEz where Xzz is the lin- 
ear dielectric susceptibility {Ez is assumed to be small). 
Therefore the polarization in the A layer is now PazE = 



Paz + ^Paz where APaz << Paz- Putting PazE in 
Eq.(36) and keeping only terms to the first power in 
APa2 leads to the equation: 

2^33^^ - Xzz [2a3 + 12033^1] + 1 = 0; T < T^f 

(37a) 

Eq. (37a) defines the static dielectric susceptibility for the 
case of inhomogeneous polarization Paz{z)- Homoge- 
neous polarization PazO leads to x^zo — 2a3 + 12a33P^^Q. 
It is obvious that the Xzz{z) dependence originates from 
the z dependence of Paz{z)- When the polarization is 
inhomogeneous, it is wrong to suppose that the suscep- 
tibility in the ferroelectric phase can be found by con- 
ventional differentiation of the free energy, i.e. Xzz — 
2a3-|-12a33P|^(z) is incorrect because Xzz obviously does 
not satisfy Eq.(37a) in general. 

Wc will now discuss 

the solutions of Eq.(37a). By introducing the dimen- 
sionless variable ^3 = 
(z3 -I- Lz/2 — 2jlAz) I a/I -I- and keeping in mind that 
2a3 + 12a33Pl(z) = 2a^{\-rk\,-%k\,sv?{i^)) j {\^k\,) 
we can rewrite Eq.(37a) in the form: 

^ + Xzz (1 + fcL - 6fcLsn2(e3)) = -2 (1 + /cL) Xt 

"S3 

(37b) 

where Xt^ — 4|a3| is the thick film susceptibility. The 
homogeneous part of Eq.(37b) is known to be a Lame 
equation |2^. The general solution of Eq.(37b) can be 
written as a sum of the fundamental solutions of the ho- 
mogeneous equation plus the particular solution of the 
inhomogeneous equation. Letting — ms and denot- 
ing the susceptibility at the boundaries as Xs we obtain 















z 


^^^^y^ y 















FIG. 1. Schematic diagram of a multilayer structure. 
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Xzziis) = XsCn(^3) dn(^3) 



X.(l^m3)^ + 2xt(l + m3)^ 
(1 + m.z)E{mj) - (1 - m.z)K(m'i) \ \ - mz 



cn(^3)dn(^3) X 



1 + m3 

1 - 7713 



i;(am(^3), 777,3) 



Sn(e3)(cn2(e3) + m2dn2(^3)) 



(1 -7773)011(^3) dn(^3) y 
+2X*71^^^ f (1 + 7^3)00(^3) dn(C3) - (1 + rna) + 277738112(^3)) . 

(1- 777,3)^ V ) 



Here 01(^3) and dn(^3) are the elliptic cosine and am- 
plitude delta functions respectively. Their forms, proper- 
ties and the relations between them are given in [p^|-^ , 
1^,^. The details of the derivation of Eq.(38) are given 
in Appendix 2. It is seen from Eq.(38) that the be- 
haviour of Xzz(C3) depends on 7773, which is the polar- 
ization ratio (see Eq. (19)). This dependence is depicted 
in Fig.7a,b,c for two values of the susceptibility at the 
interfaces: Xs = 2xt and Xs = 0. The most interest- 
ing general feature of Xzz(?3) is the appearance of peaks 
as 7773 increases. The peaks become sharper and their 
maxima tend towards the interfaces as 7773 approaches 
unity, which is when the maximum polarization of the 
layer equals the thick film polarization limit. In this limit 
Xzzi.^'i) tends to xt in the major portion of the film in- 
dependent of the Xs value, as it should in a thick film 
(ck33 — > 0, PAzm Paz)- From this we see that the so- 
lution (38) gives the correct value in the thick film limit. 

For the electric field applied perpendicular to the mul- 
tilayer surface, the susceptibility Xzz of the entire struc- 
ture can be expressed as 



-L,/2+(2j + l)Z, 
1 1 f '^^3 



-L,/2+2jU 



Xzz{z3) 



(39a) 



2K{kB^) 



Xxx 



2K{k 



Bx) 



Xxx{Cl)dCl 



(38) 



(39d) 



Since the available experimental susceptibility data for 
multilayers corresponds to the case of the electric field 
applied along the x direction, we performed calculations 
on the basis of Eq.(39d). Details of these calculations are 
given in Appendix 3. They yield the following expression 
for the multilayer dielectric susceptibility Xxx- 




This result follows from the fact that the capacitances Q 
of serially connected capacitors obey the relation 1/C = 
^ 1/Ci. When the electric field is applied parallel to x 

i 

(which will interact with Pbx only) we have the result 



FIG. 2. Variation of the polarization in a unit cell of a 
superlattice with an "a/c" domain structure. 



Af-l 



-L^/2+2(j+l}l^ 



Xxx 



7^1: 



i=o 



Xxx{zi)dzi (39b) 



-L,./2+(2j + l);, 



since here the layers are connected in parallel. The inho- 
mogeneous susceptibility Xxx{£.i) of the B layers can, of 
course, be obtained from Eq.(38) by substituting zi for 
Z3, 61 for 03, /344 for 033 and ^1 for ^3. The substitution 
of the dimensionless parameters ^3 for and ^1 for zi in 
Eq.(39a) and Eq.(39b) respectively gives 




1 

Xzz 



2K{kA,) 

1 /■ d^a 



2K{k 



Az) 



Xzziis 



(39c) 



FIG. 3. Periodic inhomogeneous polarization in a multi- 
layer structure with "a/c" domains. 
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Xxx — 



K{mi) (1 - TOi)2 Wl + mi)E{mi) - (1 - mi)K{mi) 



+ 2xt((l - mi)K{mi) - 2E{mi)) 



(40) 



where Xi~^ — 4|&i|, Jtii — kj^^. Eq.(40) gives the de- 
pendence of the susceptibihty on the temperature, the 
polarization ratio, and the layer thickness via the relation 
of these quantities (see Fig. 4). In the thick film limit as 
mi 1, Eq.(40) gives Xxx Xt (see Appendix 4), which 
can also easily be obtained from Eq.(2b) (by neglecting 
the gradient terms) as Xxx = '^^'^ b / dP'^^. The suscep- 
tibility Xxx corresponds to the ferroelectric phase of the 
thin film multilayer, i.e. T < , where is the thick- 
ness induced ferroelectric phase transition temperature 
given by Eq.(30). At this temperature the polarization is 
zero at the second order phase transition, so that in the 
limit of mi we obtain (see Appendix 4) 



Xxx ^ ixs + 2x08/37r2mi, T ^ 



(41) 



Hence we arrive at the divergence of the dielectric suscep- 
tibility for T = T^, which corresponds to Ibx = tt. This 
same result can be obtained directly by solving Eq. (37a) 
or its analog for Xxx when Paz = or Pbx — 0. In par- 
ticular, and keeping in mind that in the considered layer 
structure any function should be periodic with the pe- 
riod of the structure 21 (here we simplify things slightly 
by letting Ij^ = Ig = I @J|,|l^), we obtain: 



Xzz{z3) ^X0 + 



sin(z3 -I- - 2jlAz)iXs - Xo) 
sin(L,/2 - 2jIaz) 



(42) 



where Xo ' = -2a^{T,i " T^). 

At the transition temperature T = or where 
lz,x = and Lz^x = 2Nlz.x, one can see from Eq.(42) 
that Xzz{z) and Xxx{z) become infinitely large because 
sin(A^7r — 27rj) — 0. The divergence of the susceptibility 
at the temperature of the thickness induced phase transi- 
tion is similary obtained for Xzz and Xxx after integration 
and summation of Eqs.(39a) and (39b). Therefore the di- 
vergence of the susceptibility Xzz or Xxx at the transition 



temperature T = oi T 



is the characteristic fea- 



ture of ferroelectric thin film multilayers. In Fig. 8a we 
display the temperature dependence of the susceptibil- 
ity in the ferroelectric region T < for the thickness 
Ib — 21bo, (recall that Ibo is the thickness where van- 
ishes) which corresponds to /T^x = 0-75 (see Eq.(30)) 
and Xs = 0- Since Ib — Ibo gives = this ratio 
will be smaller than 0.75 at Ibo < < 2^bo and larger 
than 0.75 at Ib > 21bo- K appears possible to write the 
temperature dependence given by Eq. (40) in the form of 
a Curie- Weiss (C-W) law at T ^ T| (mi -^0). We 
obtain from Eq.(41) the following approximate forms of 

Xxx 

(r)for three different values of x^: 



Xxx{T) 



1 



cb 



r2nB rpB 



T"-T 



^0 ^cl 



rpB 



Xs = (43) 



Xxx{T) 



4y ^ 

,2 

f- r 



T^l - T) 



Xs = const 7^ 

Xxx(T) w 
Xs = axt 



2 + a 



a 

" 4ai {TB-T) 



(44) 



(45) 



Details of these calculations are given in Appendix 5. 
Note that the C-W constant for Xs = is very close to 
that of the thick film value, Ct = (4a^)^^. It is seen from 
Fig. 8a that Eq.(43) fits surprisingly well the temperature 
dependence given by Eq.(40) (see solid curve in Fig. 8a) 
not only in the vicinity of T = but also for the entire 
temperature region. Eq.(45) also gives a good fit after 
renormalization of the Xi value in Fig. 8a. Eqs.(43) and 
(45) show that the susceptibility dependence on the film 
thickness is mainly defined by the thickness dependence 
of . We point out that for xs = const ^ 0, Cb is also 
thickness dependent (see Eq.(44)). 

Well into the paraelectric phase at T > or T > T^f 
where there is zero spontaneous polarization, there will 
be some small homogeneous polarization induced by the 
external electric field. In this case the free energy can be 
expanded in a power series of the polarization, and the 
static dielectric susceptibility can be found, as usual, by: 



Xz 



^2$ 

ap2 



(46) 



Eq.(46) leads to conventional C-W laws for the suscepti- 
bilities in the paraelectric phase of thin films which sup- 
plement the susceptibility in the ferroelectric phases of 
the "A" and "B" layers 



Xzz 
Xxx 



t>t:, 



T>Ti 



(47) 
(48) 



which diverge at the transition temperatures or 
which can be considered as the Curie temperatures of 
the thickness induced ferroelectric phase transition with 
the characteristic dependence on the film thickness (see 
Fig.5). Note that for thick films, when ~f and 
(see Fig.5), the coefficients Cq and — > 

Uq , but in general cq and ao are different. Unfortunately 
it is cumbersome to compare the values of cb in Eq.(43) 
with those of Cq in Eq.(48). For the sake of illustration 
we depict in Fig. 8b the Xxx dependence for the entire 
temperature region in supposing that cb — 1/{2cq) and 
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Cb = The divergence at T = is the most 

significant feature here. 

It will be shown later, by comparison with the exper- 
imental data (7), that the giant dielectric response ob- 
served in some multilayers originates from the suscepti- 
bility anomaly related to thickness induced phase tran- 
sition considered above. Of course the internal fields in- 
duced by misfit dislocations, growth imperfections and 
impurities will smear out the dielectric response of the 
multilayers and thus reduce infinity to a more realistic 
value. 

We will now discuss the possibility of calculating the 
thin film susceptibility as the second derivative of the 
free energy instead of solving the Lame equation. To do 
this we assume that in the ferroelectric phase the role of 
the second term in Eq.(37) will increase with increasing 
polarization, and so it should be possible to neglect the 
contribution of the second derivative. In this case the sus- 



ceptibility becomes Xz 



2a3 -I- Integration 



of this expression over the layer thickness (with Paz{z) 
given by Eq.(23)) and then summing over the layers leads 
to the following expression 



= 2a3 

Xzz 



6 



1 



EjkAz) 
K{kAz) 



(49) 



where Uaz is given by Eq.(19). 

Equation (49) gives the correct expression for the sus- 
ceptibility in the thick film paraelectric phase {kAz = 0, 
X^} = 2a3, T > T^) and in the thick film ferroelec- 
tric phase {kAz = 1, Xji = — 4a3, T < T^). However, 
problems arise with this approximation for intermediate 
values of kAz- For kAz < k\^ = 0.674 we obtain a nega- 
tive expression for Xzz (including the region of thin film 
paraelectric phase T > T^) and a positive expression 
when 1 > kAz > k\^ with Xzz diverging at k^z — 0.674 
i^Az sends the expression (49) for x7z to zero). It is ob- 
vious that for kAz < k\^ the susceptibility must be calcu- 
lated from the differential equation Eq. (37) instead of the 
approximation derived Eq.(49) since the latter leads to 
the physically unreasonable negative susceptibility. The 
validity of Eq.(49) in the region fc^^ < kAz < 1 when 
1/Xzz > can be checked by a 1/xzz calculation on 
the basis of the exact solution of Eq. (37) . Thus we see 
that Eq.(49) will be valid only in the thick film limit of 
kAz = 1- Therefore, for the majority of thin films and 
multilayers, the dielectric susceptibility will be defined 
by Eq. (37) whereas the normally used free energy second 
derivative (see, e.g. jl|]) is valid in the ferroelectric phase 
only for thick films. We also point out that Eq.(37) will 
be applicable to bulk ferroelectrics with inhomogeneous 
polarization, i.e. when there is a strong contribution of 
the polarization gradient since the mean field approxima- 
tion is no longer valid. 



III. DYNAMIC PROPERTIES 

A. Polarization 

To consider the dynamic properties of a periodic spa- 
tial structure, we should add the time derivatives of the 
polarization to equations (12a) and (12b]. This can be 
done by the procedure suggested in [P5| , p6| . This proce- 
dure permits the investigation of low frequency dynamics 
in ferroelectrics. Its essence is to expand the polarization 
in powers of frequency or the time derivative operator 
d/dt. The first order term gives, as usual, the decay of 
the polarization, while the second order term provides an 
"oscillatory response" and contains the mass coefficient 
^. 

The equations of motion can be written as 



d'^PAz , OPaz 
t^A^rp^ 1- lA- 



fJ^B- 



d^PB. 



dt 



SFa 
SPaz 



= 



dP, 



7b- 



Bx 



dt 



SFe 

SPbcc 



= 



(50a) 



(50b) 



The free energy variation SF/ SP in the Euler-Lagrange 
Eqs.(12a) and (12b), yields the equations for Paz and 

Pbx 



^ ^'"^ ZD 0.2 0.4 0. 




0.0 0,2 



FIG. 4. Dimensionless thickness of an "A" or a "B" layer 
versus the polarization ratio. 




FIG. 5. Temperature of the thickness induced ferroelectric 
phase transition as a function of layer thickness. 
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dPAz ^ p 

at 



dP, 



7s- 



Bx 



dt 



Since Eqs. (51a) and (51b) are essentially the same 
we will consider the solution for Paz keeping in mind 
that Pbx can be written, as usual, by substitution of the 
relevant coefficients and constants in Eq.(51). 

We will look for solutions of Eq.(51a) in the form of sta- 
tionary nonlinear waves in the conventional self-similar 
form 

Paz ^PAz{z-vt)=P,{kz-ujt), uj^kv (52) 

where uj, k and v are the frequency, the wave vector, and 
the velocity respectively and both vectors are directed 
along the z axis. Since the wave propagates along the 
normal to the multilayer structure (fc and v are parallel to 
z) we have suppressed the index z in k and v in Eq.(52). 
Using the relations 



9*2 



17' al2 = ^ 172' = 1P2 (^'^) 



we obtain from (51a) 

. d^PAz 



iHAV -033)- 



de 



lAV- 



dPAz 
d^ 



(54) 



Since we are interested in the dispersion law for non- 
linear waves (which aX t — Q gives us our static periodic 
structure (23)) we shall neglect the decay term. In doing 
so (7A — 0) Eq.(54) coincides with Eq.(12a) but with 
ck33 — > a33 — ^jlav^ , so that the solution of Eq.(54) can be 
written in the form of Eq.(23) : 



Paz{z - vt) ^ Paz 



iSn 



zz- vt + Lz/2 - 2jlAz 



^Az 



, ^Az 



(55) 



where v — \/ —cij,/ {a^-^ — fj,Av'^)v and the renormalized 
values of Z3, L2, and Iaz can easily be obtained from (22) 
by substituting ^7-03/(033 - fiAv'^) for ^^3/033. 

The solution of Eq.(51b) can be similarly obtained and 
written in the form: 



PbxIz - Vt) = PBxmS.Tl 



Zi -Vt + La;/2-2jlBa 



, kBx 



(56) 



y/-bi/{(344 - Hbv'^)v and zT, L^, Iaz 
Eq.(29) after the substitution of 

(55) and (56) 



with V 
follow from 

a/-6i/(/344 - y-Bv'^) for ^J-bi/j3ii. Eqs. 
give the time-dependent inhomogeneous polarization of 
a multilayer structure. 



blPBx + 2611PI2; ^ A 



O^Pbx 
^ dz"^ 



: 



= 



(51a) 



(51b) 



B. Dispersion law 



To find the dispersion law of the stationary nonlin- 
ear waves formed from Eqs. (55) and (56) in the A 
and B layers respectively we first recall that the wave 
number k = ^njiXAz + ^Bx)^ where \az and \bx are 
the wavelengths of the Paz and Pbx nonlinear waves 
respectively. They are actually determined by Eqs. (25) 
and (28b) with respect to the substitution a33(or ji^^) — > 
a33(or /?44) - lJ'A.,Bv\ g. Hence 



' ^33 


- HAV'^ 


-as 




- fJ-BV^ 



1 + k\^K{kAz) 



^Bx ~ - 1/ , 

It then follows that 



1 + kl^K{kBx) 



(57) 



3.5 4 5 6 7 /„. 




0.2 0.4 0.6 0.8 1 kBz,x 
0.20.4 0.6 0.8 1 ^B^' 



FIG. 6. Form of the function tp{k) in Eq.(32) which deter- 
mines the criterion for the multilayer domain structure. 
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27r 



TT 



Expression (58) implicitly determines the dispersion 
law, since v = uj/k, and so Eq.(58) can give the uj{k) 
dependence. Let us simplify Eq. (58). It follows from 
Eqs. (25) and (28b) that 



Ia 



2 

5=2 



Q33 
-03 

/344 
~1 



k\^K{kAz) 



1 + kl^KikB.) 



Thus we can rewrite Eq.(58) as: 

1 



k = 2tt- 



(59) 



(60a) 



Xzz(^3,m3) 

















1113=0.2 




















1113=0.3 


























m3=0.5 













































0.2 0.4 0.6 0.8 



2K(m.,) 



(58) 



or 



"33 



MS 2 

P44 



27r 



(60b) 



To facilitate the analysis of the dispersion law given by 
Eqs. (60) we will first assume that rj = tja = Vb, where 
fJ-A/aas = VA, (J-b/Pm = Vb and Ia ^ h ^ I- In this 
case we have from Eq.(60): 



k^ - — 



(61) 



It is seen from Eq.(61) that the dispersion law has the 
familiar long wavelength form. It also exhibits the pe- 
culiar dependence of the nonlinear wave frequency on its 
amplitude via relation between I and kz x (Eq.(59)), see 
e.g. H. 

We schematically plot the dispersion law in Fig. 9a. 
First of all since cj^ > we have some critical value 
k = kc = ±7r// for which w = . 

On the other hand for any specific value k = kg there 
is a critical thickness Ic = ±7r/fco at which lu — and 
uj for I > Ic only. The thickness dependence of the 
frequency has the form: 



Iko 



(62) 



Ic gives u! 



which in the vicinity of / = 
±y/l/r]koy/2/lcy/l — Ic- Thus the frequency increases as 
the square root of I . This dependence is represented in 
Fig.9b. 

For the general anisotropic case, rjA 7^ lys and toil) 
and uj{k) will be represented by more complicated ex- 
pressions. They can be derived by squaring both sides of 
Eq.(60b). We obtain ior U ^ h ^ I- 



Xzz(?3,m3) 



2Xt 



5Zt 
Xzz(ii,m3) 



3Xt 



0.2 0.4 0.6 0.8 



















013=0.7 












^ — m3 


1113=0.9 
= 1-10"^ 



















2K(m3) 







m3=0.2 




















113=0.5 














m3 


=1-10"'' 

















1 2K(m3) 



FIG. 7. Inhomogeneous susceptibility dependence on the 
coordinate ^3/(2A'(m3)) — zi/Iaz for different ms values for 
Xs = 2xt (a and b) and for xs = (c). 



12 



We have dropped the minus sign before the square root 
in Eq.(63) because ajf 2 > 0- K is seen that when — 
TT^/P then u!i^2 = 0, similar to the isotropic case. Thus 
this noteworthy aspect of the dispersion law is conserved 
for the general anisotropic case. We plot the general 
dependence in Fig. 10. Although the qualitative features 
of the uj(k) law for the isotropic case are preserved in the 
anisotropic case, there are peculiarities at large k in the 
anisotropic case: we determine uj ~ \/k rather than linear 
dependence in the isotropic case (compare Figs. 10a and 
9a) and there is a different thickness dependence, with a 
small decrease of oj{l) at large I (compare Figs. 10b and 
9b). 



1 

(a) ^^^^ 









0.2 0.4 0.€ 0.3 1 

T 
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1 
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1 
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0.2 0.4 0.6 0.8 1 

T 



FIG. 8. Temperature dependence of the dielectric suscep- 
tibility in units of xi = l/{al^T^^), for xs = (a) at T < 
using Eq.(40) (solid line) and Curie- Weiss law (dashed line) 
for Ib ~ '21 BO and (b) in the entire temperature range using 
Eq.(40) (T < rj) and Eq.(48) (T > TJ*) for different vab 
ues and for Ib ~ 21bo (solid line: — 2cs; dashed line: 

1/C|f = cs/2) and Ib = SIbo (dot-dashed line : 1/cq = 2cb\ 
dotted line: 1/C(f = Cfl/2). 



(a: 



CO 




7i/k,pl, 



FIG. 9. Nonlinear wave dispersion law for a multilayer in 
the isotropic case as a function of (a) the wave vector and (b) 
the film thickness. The dashed line is a; = kri~^^^ . 



IV. DISCUSSION. GIANT DIELECTRIC 
RESPONSE OF FERROELECTRIC THIN FILM 
MULTILAYERS. 

Let us begin with the comparison of the calculated and 
the observed dielectric susceptibility. The susceptibil- 
ity was recently measured [0 in superlattices consisting 
of ferroelectric PbTiOa and paraelectric Fbi-^^LaajTiOa 
(x = 0.28 - PLT) grown on (lOO)-oriented SrTiOa single 
crystal substrates. Three superlattices of PT/PLT with 
modulation wavelengths of 100 A° (sample S-40), 400 
A° (sample S-10), 2000 A° (sample S-2) were studied. In 
each superlattice the FT and PLT layers were of equal 
thickness, and the total thickness was 4000 A°. The au- 
thors observed a Debye-like frequency dispersion of the 
real and imaginary parts of the dielectric susceptibility 
Sxx for the samples S-10 and S-2 in which the value of the 
real part of the dielectric susceptibility at low frequency 
approached 420000 and 350000 respectively at T w 50°C. 
A significant increase in the susceptibility with tempera- 
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ture was observed in sample S-10. We present this tem- 
perature dependence by the full circles in Fig. 11, where 
the solid line represents the C-W law given by Eq.(43) 
with CB = 7 • lO^K, and T^f (S-10) = 533K. The inde- 
pendence of cs on layer thickness (see Eq.(43)) makes it 
possible to obtain the transition temperature for speci- 
men S-2: T^f (S-2) = 575K. Taking the ratio of these 
values and using Eq.(30) for S-10 and S-2 leads to 
the critical thickness I bo = 55 A°. This value is larger 
than the layer thickness {Ib = 50yl°) in S-40 specimen, 
which is why no ferroelectric phase transition was ob- 
served for this multilayer. Thus the value reported for 
S-40, e'^^ ~ 750 at T « 50°C, represents the contribu- 
tion from the paraelectric phases of both the PT and PLT 
layers. This value is several orders of magnitude smaller 
than the values measured in the ferroelectric phase for S- 
10 and S-2. The contribution from the PT-layers in the 
paraelectric phase of the S-40 sample, where Td ~ (see 
Eq.(30) at /s ~ Ibo) can be written as Sxx = "^t^cbp/T, 
where CBp — Cq in accordance with Eq.(48). Keeping 
in mind that e'^^iPT) < 750, wc can conclude that the 
C-W constant in Eq.(48) csp < 2 • lO". Therefore the 
dielectric susceptibility in the paraelectric phase of the 
thin film is much smaller than in the ferroelectric phase. 
This is contrary to what is observed in bulk ferroelectrics. 
On the other hand this gives support to the supposition 
that the contribution of the permittivity from the PLT 
paraelectric phase to the giant value observed in S-10 
and S-2 is negligibly small in comparison with the con- 
tribution from the PT layers. Estimation of the transi- 
tion temperature for a thick layer, T^, with the help of 
Eq.(30) and the values determined above for leads 
to — 580K. This value is smaller than that for bulk 
PbTiOa : Tjj = 763K. This would seem to favor com- 
pressive strain in the layer, i.e. < in Eq.(9b) be- 
cause Qfi + Qfa > for PbTiOg. Note, that the 
and cb values were obtained assuming the temperature 
dependence of the low frequency (1 kHz) dielectric per- 
mittivity is close to the static one. Calculation of the 
static susceptibility via the maximum value of e^x mea- 
sured at T = 50°C has shown that the low frequency 
e^^ was in fact close to the static value. 
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FIG. 10. Nonlinear wave dispersion law for a multilayer in 
the anisotropic case at rjBx = 0.2, tjbz = 0.25 (1); rjAx ~ 0.2, 
rjAz = 2 (2) as a function of (a) the wave vector and (b) the 
film thickness. 



V. CONCLUSION 

The model used for the calculations was a multilayer 
consisting of alternating layers of two different ferroelec- 
tric materials for which the interfaces between the layers 
were abrupt and the polarization at these interfaces was 
taken to be zero. The interaction between the polariza- 
tion in different layers was neglected because of the rapid 
decay of this interaction with distance jl^ and because 
of the zero polarization at the interfaces. The produc- 
tion of sharp interfaces in multilayer structures is now 
possible using the most recent oxide thin film deposition 
techniques as has been demonstrated in . Our indepen- 
dent layer model produces results that are in agreement 
with published results |^,^,|l^ both for multilayers and 
for single thin films. In addition, this model should also 
be valid for bulk materials in which polarization gradi- 
ents are present and thus rendering calculations in the 
mean field approximation meaningless. 

We have shown that the inhomogeneous polarization 
in the multilayers is at the origin of the inhomogeneous 
static dielectric susceptibility. We have determined the 
differential equation for the susceptibility and, from its 
solution, we show that the conventional way of calculat- 
ing the susceptibility as the second derivative of the free 
energy jl^ is approximately valid only for thick films 
with small polarization inhomogeneity. The results re- 
lated to the consideration of renormalization of the free 
energy coefficients by homogeneous stresses in the film, 
to the criteria of "a/c", "c/c" or "a/a" domains struc- 
ture as well as to thickness induced ferroelectric phase 
transition can be applied both to multilayers and, after 
some transformations, to single thin films on substrates. 
Several results are valid mainly for ferroelectric thin film 
multilayers. In particular, the results obtained for the 
nonlinear polarization waves and for the dispersion law 
exhibiting a critical wave vector or, equivalently, a criti- 
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cal thickness, reflect the characteristic peculiarities of the 
dynamic properties of thin film multilayers. 



VI. APPENDICES 
A. Appendix 1 



e'xx, xlO= 



T, °C 



40 60 80 100 120 140 



FIG. 11. Temperature dependence of the dielectric permit- 
tivity of a PT-PLT multilayer. Solid line: theory; Full circles: 
experiment [7]. 

The calculations of the inhomogeneous susceptibilities 
in the form of Eq.(42) can be valid only for multilayers. 
Therefore the calculated divergence of x at the thickness 
induced ferroelectric phase transition temperature may 
be the exclusive feature of ferroelectric thin film multi- 
layers. Because of the technical difficulties in producing 
and studying ferroelectric thin film multilayers are more 
involved than for those of single films or bulk materials, 
experimental results are still somewhat restricted to just 
a few different series of superlattices. The theoretical cal- 
culations of multilayer structures performed in this work 
are quite complex even for the relatively simple indepen- 
dent layer model. This model corresponds fairly closely 
to that of ||l^ where, for the case of weak coupling at 
the interfaces, they also find the most interesting phys- 
ical phenomenon in the layers: a thickness induced fer- 
roelectric phase transition. Our results are in agreement 
with the numerical solutions in Ref. [|l2|. We are now 
extending our calculations to obtain the static and dy- 
namic dielectric susceptibilities over a wide temperature 
region, including the response of the ferroelectric phase 
for different orientations of an external electric field. A 
more general model should take into account the inter- 
faces where gradients in the mechanical stress can appear. 
Thus the calculations will increase in complexity since 
the equilibrium stresses must then be calculated on the 
basis of the entire Euler-Lagrange equation (7b) rather 
than by the simple differentiation of the free energy with 
respect to stress that we have used above. Nevertheless 
even with the increased computational complexity, fur- 
ther theoretical, as well as experimental, investigation of 
ferroelectric multilayers, including those with finite in- 
terfaces and with more than two layers in a modulation 
period, would seem to be extremely desirable both for 
basic science and potential applications. 



Since the free energy density fsiPsz) and fBiPsx) as 
well as Pbz and Pbx can be obtained from one another 
by substitution of their parameters (see Eq.(5a) A ^ B 
and (5b), Eq.(23) with (28a)) we shall perform a detailed 
calculation of Fb{Pbz) and than obtain Fb{Pbx) by this 
substitution. 

To calculate Fb{Pbz) we can take into account the first 
integral for Pbz (the step going from Eq.(14) to Eq.(15) 
with A ^ B). In view of this first integral we can sub- 
stitute the derivative for powers of Pbz in Eq.(5a) (again 
with A ^ B) and obtain 



fB{PBz) = 2(3: 



'33 



dPsz 
dz 



CBz 



cbz = bsP'i 



3^Bzm ^33^B 

We determine the derivative to be: 

dPp. ' 



(Al.la) 



(Al.lb) 



dz 



P, 



1 



Bzrn 1 



/333 1 + 



cn{u^) dii{uz) (Al.2) 



z + L/2-i2j + l)lB 

>33 V^TWz 



where cn(u) and dn(M) are elliptic functions (see e.g. 
]l9[). With respect to (Al) and (A2) the free energy 
Fb{Pbz) can be written as 

^ -L/2+2(] + 1)Ib . 

Fb{Pbz)^\Y. J dzlcBz- (A1.3) 



-P, 



2k 



'3 „„2/, N ,„2 



Bzm ^ _^ p 



cn {uz)dn (u^) 



Bz 



h 



CBz -Pi 



2 ^"3 

Bz 



2h 1 



y/TT^jBz 



X / cn^(u) dn^(u) > du. 



Since Ibz = 2^1 + k%^K{kBz) (see Eq.(25) with A 
B) the last integral in (A3) has the form 

2K{kBz) 

1= J cn^{u)dn{u)du (Al.4) 



We will now calculate the integral (A4) (see p^) 
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2K{kB^) 



1 = 2 J (1 - sn^u){l - kBzSn^u)du = 



= 2{K{kB.)-{l + k%Jh+k%j2}, 



(A1.5) 



where 



/i = J sn^udu, 




I2 = I sn^u du 







^cn(M) dn(u) sn(u) 1^^'=^^^ + 



2(i + fcL) 



K{kB.) 



/ sn^udM- -l-ii:(fcB^). (A1.6) 



Since cn{K) = sn(0) = 0, we obtain 



T 1 K(k ^-^2a + 4J 



(A1.7) 



We calculate /i 



/i = {^i - £;(am(«), fcBz)} 1^^'"^^ = 

= T^{K{kBz)~E{7T/2,kBz)} = 



1 

''^Bz 



7r/2 



(Al. 



where ^^(fcBz) = J wl — fc^^sin 0d0 is complete ellip- 
* 

tic integral of the second kind. 

Taking into consideration (A1.5), (A1.7), and (A1.8), 
the integral (A1.4) can be rewritten as 



2K{kBz) 



^l±^{K{kBz)-E{kBz) 



(A1.9) 

Substitution of (Al.9) and (Al.la) into (A1.3) gives 



FBiPBz)= 



'■B^Bzm 
I A + Ib 



h + bssPBzm 



{l + kl^)K{kBz) 



2K{kBz) - 



1+3 

^Bz 



'-{K{kBz) - E{kBz)) 



(ALIO) 



To express (Al.lO) in terms of the parameter kBz and 
the coefficients &3, 633 we first substitute Psztti/^BzO — 
2^Bz/(l + ^Bz) (see Eq.(19)) and then we take into ac- 
-63/(2&33). Thus we have 



count that -Pj^o 



Fb{Pe 



1 



X % + l-(l 



Ia + Ib b33 {I + kEz)"^ 3 

E{kBz) \ 



2 

- X 



klz 



K{ki 



(Al.ll) 



and Fb{Pbx) can be obtained from (Al.ll) by substitut- 
ing kBx for fcsz, 61 and 611 for 63 and 633 respectively. 
This gives the inequality (32) and the function (33). 



B. Appendix 2 

Let's consider the solution of Eq.(37b) for the dielectric 
susceptibility: 



+ (1 + TO - 6msn^(C,TO)) xiO = -2(1 + TO)xt, 



with the boundary conditions 

x(e - 0) - xii - 2i^(TO)) = xs 



(A2.1) 



(A2.2) 



Here for the sake of simplicity we have omitted all sub- 
scripts and have denoted m = ^bx- accordance with 
the theory of linear differential equations we have to solve 
the homogeneous part of equation (A2.1), which is a par- 
ticular case of the Lame equation |^,^ . A solution of 
the homogeneous equation is the following [M: 



XiiO = Ci cn(^, m) dn(^, m), 



(A2.3) 



A second solution can be easily obtained by varying the 
constant Ci and it has the form: 



X2{i) = C2 cn(^,m) dn(^,m) / 

Jo 



dc 



cn2(C,m) dn (C,^) 

(A2.4) 



The particular solution of the inhomogeneous equation 
(A2.1) can be obtained by varying the constants Ci and 
C2 in the sum of the solutions (A2.3) and (A2.4). This 
procedure leads to: 



XiiO ^ -2(1 + m)xt cn(^, m) dn(^, m) x 
sn(C, m) dC 
cn2(^, m) dn^(C, m) 



(A2.5) 



After substitution of the variable ( hy (p = am(C), this 
integral has the form /i = Jo"^^^\^ ~ sin((/3)^)~^(l — 
msm{ip)'^)~i . The integral /i can be put in standard 
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form (see ||T^), and its substitution into Eq.(A2.4) gives 
us the direct form of the second fundamental solution of 
the equation (A2.1): 

/ 1 + m \ 
X2{0= U - -j -B(am(^),m) cn(,^, to) dn(C, m) + 

\ 1 — TO / 

^ sn(^,m) (cn^(C,TO) + TO^dn^(^,TO)) 

1 + TO 

The integral in the expression (A2.5) can be reduced, 
by the substitution t = cn(C,TO), to that which can 
be found in various sources pH.Then Eq.(A2.5) can be 
rewritten as 

X%{i) = ^Xt- ^2 (1 + "^) cn(^,TO)dn(^,TO) + 



2K{m) 



(1 — mY 



+2TOSn^(^, to) — (1 + to) 



(A2.7) 



Finally the complete solution of Eq.(A2.1) has the fol- 
lowing form: 



x(0 -Cixi(e) + C2X2(0 + x.(0 



(A2., 



where the constants Ci , C'2 are determined from the 
boundary conditions (A2.2). When using the proper- 
ties of elliptic functions , the constants Ci , C2 can be 
found to have the form: 

Ci - xs, (A2.9) 

(l-TO)2xs+2xt(l+TO)2 f 1 



(1 + m)E{m) — (1 — m)K{m) \1 — to 



Substitution of (A2.3), (A2.6), (A2.7) and (A2.9) in 
(A2.8) gives us the expression (38). 



C. Appendix 3 

We will calculate the susceptibility from Eq.(39d) 
where, for the sake of simplicity, we have omitted all 
subscripts and have denoted to = fc^^: 



2K(m) 



X{m) 



2K{m) 



using Eq.(38) for x(Ci'^)- It is convenient to divide the 
integral in the expression (A3.1) into the following three 
parts: 



2K(m) 



2K{m) 



J Xi(C)de= J cn(Odn(e)de = sn(e) 



— ^' 



(A3. 2a) 



X2{m 



2K(m) 

l^LJH I sn(n(l + TO-2TOSn2(n)d^ 

1 — TO J 


1+"^ li:\A fc\\'^K{rn) 1 + TO 

-Y-77^™(^)dn(C) lo 



1 — TO 



(A3. 2b) 



2K{m) 



((1 + to) (cn(^) dn(e) - 1) + 2TOSn2(^)) 



2((1 - m)K{m) - 2E{m)). 



(A3.2c) 



Grouping equations (A3. 2a), (A3. 2b) and (A3. 2c), one 
can obtain Eq.(A3.3) which coincides with Eq.(40): 

_ 1 1 -I- TO 
^^""^ " K{m) (1-to)2 

^ / {1-m f XS + 2X1(1 + m f 
^1(1 + m)E{m) - (1 - m)K{m) 

+2xt{{l~m)K{m)-2E{m))^. (A3.3) 

D. Appendix 4 

It is necessary to consider two extreme cases of (A3. 3). 

E. Case 1: m ^ 1 (thick film limit). 

When we expand the elliptic functions as a series about 
fi = [1 — to) ^ we obtain: 



2 

IT f ^ 1 3 2 



and using the Lcgcndrc relation |T9| 



(A4.1) 



(A3.1) E{m)K(l -m)+ E{1 - m)K{m) - K{m)K(l - to) 



we find that: 



and 



(1 -I- m)E{m) - (1 - m)K{m) = 



2E{m) - (1 - m)K{m) = 



(A4.2a) 



(A4.2b) 
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After substituting Eqs.(A4.2a) and (A4.2b) into the ex- 
pression (A3. 3) and keeping in mind that K(m) ^ oo as 
m ^ 1 we find the following limiting expression: 

X(m ^ 1) = 

= {xt + [xs + j^) (l - - Xt 

(A4.3) 



From this temperature dependence of mi we can find 
from Eq.(A4.4) the susceptibility Xxx as a function of 
temperature. In order to obtain the temperature de- 
pendence of the susceptibility we have to choose a spe- 
cific susceptibility Xs{T) at the interface. We considered 
two forms for Xs on temperature and obtained the corre- 
sponding formulae for Xxx{T) : 

Xs = const 



F. Case 2: m — » (paraelectric thin film limit). 

For m ^ the thickness of the layer tends to the 
critical value Izx = tt at which the transition to the 
paraelectric phase occurs. Since at to ^ 0, K{m) 
(7r/2)(l + m/4) and E{m) (7r/2)(l - to/4) (which can 
be seen from the integrals of Eqs. (24) and (A1.8)), we 
obtain from Eq.(A3.3): 



X(m ^ 0) ^ {xs + 2xt) 
which coincides with Eq.(41). 



377^ TO 



(A4.4) 



Xx.(T)«(4x.§i + ^ 



Xs = aXt 



1 1 



Xxx{T) 



(A5.5a) 



(A5.5b) 



Under the assumption Xs = 0, Eqs.(A5.5a) and (A5.5b) 
lead us to Eq.(43). 



G. Appendix 5 

Let us find the approximate dependence of the suscep- 
tibility (A3. 3) on temperature. Only the dimensionless 
thickness Ibx depends on temperature in our theory (see 
Eq.(29)): 
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(A5.1) 



In the vicinity of the thickness induced transition tem- 
perature T = Tci{Ibx = "■) the expression (A5.1) can 
rewritten as: 



Ibx 



'BO" 



(A5.2) 



where Ibo and Td are given by Eq.(30). 

On the other hand it follows from the expansion of 
Eq.(25) at TOi ^ that: 



Ibx = 2Vl + TOiif(TOi) TT [1 + -nil 



4 fhx , 
mi^ -[ 1 

3 \ TT 



Substitution of Eq.(A5.2) into Eq.(A5.3) yields 
2 1% - T 

'bo -^cx 



(A5.3) 



(A5.4) 
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